Thermodynamics of the bilinear-biquadratic spin one Heisenberg chain 
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The magnetic susceptibility and specific heat of the one-dimensional 5 = 1 bilinear-biquadratic 
Heisenberg model are calculated using the transfer matrix renormalization group. By comparing 
the results with the experimental data of LiVGe2 0e measured by Millet et al. (Phys. Rev. Lett. 
83, 4176 (1999)), we find that the susceptibility data of this material, after subtracting the impurity 
contribution, can be quantitatively explained with this model. The biquadratic exchange interaction 
in this material is found to be ferromagnetic, i.e. with a positive coupling constant. 



PACS Numbers: 75.10.Jm, 75.40.Mg 

The quantum spin chains have been the subject of 
many theoretical and experimental studied since the con- 
jecture made by Haldane 0] that the antiferromagnetic 
Heisenberg model has a finite excitation gap for integer 
spins. The model which has been intensively used to in- 
vestigate the physics behind the Haldane's conjecture is 
the isotropic spin one Heisenberg Hamiltonian with both 
bilinear and biquadratic spin interactions: 



H = JV[S,.S 1+1 + 7 (S,-S 1+1 ) 2 ] 



(1) 



For most of the existing quasi-one-dimensional (ID) 5 = 
1 materials, the biquadratic term is very small compared 
with the bilinear term as well as the uniaxial anisotropy. 
This model was therefore generally thought to be of pure 
theoretical interest. However, recently Millet et al. 
found that the magnetic susceptibility of a new quasi- ID 
5 = 1 system, the vanadium oxide LiVGe20g, shows a 
few interesting features which are absent in other 5=1 
materials. They argued that both the interchain coupling 
and the uniaxial anisotropy are too small to create these 
features and suggested that the biquadratic term plays 
an important role in this material. 

In this paper, we present a theoretical study for the 
thermodynamics of the bilinear-biquadratic spin chain 
(Q) with J > 0. We have calculated the magnetic suscep- 
tibility and specific heat of this model using the transfer 
matrix renormalization group (TMRG) method |^-|(| . By 
comparing with the experimental data of LiVGe2 06, we 
find that the measured susceptibility, after subtracting 
the impurity contribution, can be quantitatively fitted 
by the numerical result with 7 = 1/6. This shows that 
the spin dynamics of LiVGe 2 C>6 can indeed be described 
by the Hamiltonian ([[]), in agreement with Millet et al. 

However, the value of 7 needed for fitting the exper- 
imental data is different from that suggested by Millet et 
al. @. 

Let us first consider the properties of the ground state. 
It is known that when 7 = — 1 and 1, the model ([!]) can 
be solved rigorously by the Bethe Ansatz Between 
these two soluble points, the system is in the Haldane 
phase. In this phase, the ground state is a non-magnetic 



singlet with a finite energy gap in excitations. In par- 
ticular, when 7 is between —1 and ji C ~ 0.41, the low 
energy physics of this model can be understood from the 
valence bond solid (VBS) model proposed by Affleck et 
al. P]. In this model, each site on the chain is oc- 
cupied by two 5=1/2 spins and the ground state is 
formed by the bonding of two 5=1/2 spins from adja- 
cent sites. These singlet bonds must be broken in order 
to excite the system and this leads to a non-zero energy 
gap in the low-lying spectrum. This picture has been con- 
firmed experimentally jl^,[ll| as well as numerically |L2"| . 
At 7i C , the ground state undergoes a commensurate- 
incommensurate transition and the critical exponent for 
the magnetization changes from 1/2 below ji C to 1/4 
above 7; c |T^,[l4|]. Between j ic and 1, the system is in 
the incommensurate phase and the incommensurate peak 
in the spin form factor S(q) of the ground state moves 
continuously from 7r to 2tt/3 as 7 increases from ji C to 1 
Jl5| , ^6[ . Above 7 = 1, the true nature of the ground state 
is still controversial p^l7[|. Some works |T5|Jlt^ ] suggest 
that it might be in a trimerized phase. When 7 < — 1, 
the ground state is doubly degenerate and dimerized. 

The TMRG is a finite temperature extension of the 
powerful density matrix renormalization group method 
|]l8f . A detailed introduction to this method can be found 
in references The TMRG method handles directly 

infinite spin chains and thus there is no finite system size 
effects. To calculate the spin susceptibility, we first eval- 
uate the magnetization M of the system with a small 
external field B, and then from the ratio M/B we obtain 
the value of the susceptibility. The specific heat is evalu- 
ated from the numerical derivative of the internal energy 
with respect to temperature. At low temperatures, since 
the specific heat is very small, the relative error of the 
specific heat may become quite large. In most of our 
calculations 100 states are retained. 

Figure [l] shows the zero-field spin susceptibility x(^) 
normalized by the its peak value x pea k as a function 
of the normalized temperature T/T£ cak for a set of 7, 
where Tp 0ak is the temperature of x p0 ak. Above Tp eak , 
x(T)/x pea k behaves similarly for all the curves shown in 
the figure. When 7 is positive, x(^) drops quickly below 
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Tp Cak . This is because the energy gap in this parameter 
regime is very large. As 7 becomes negative, xC^O J us t 
below Tp eak tends to become flatter. At 7 = — 1, there 
is no gap in the excitation spectrum, x(T) shows a small 
positive curvature at low temperatures, as in the S = 1/2 
Heisenberg chain. 
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FIG. 1. The normalized spin susceptibility x/Xpcak as a 
function of T/Tpeak at zero field. The experimental data of 
LiVGe2 06 obtained by Millet et al M (empty circles) are also 
shown for comparison. The inset shows the 7 dependence 
of the peak susceptibility Xpeafc (empty circles) and the peak 
temperature Tp Cak (filled circles). J is set to unit. 

The inset of Fig ure ^ shows the 7 dependence of Xpcak 
and T p s cak . The increase of Xpcak with 7 indicates that 
the susceptibility becomes larger when 7 moves from the 
dimerized phase to the Haldane phase. This is consistent 
with the picture that in the dimerized phase the spin 
is frozen by forming rather rigid spin singlet, while in 
the Haldane phase the spin is relatively free above the 
Haldane gap. The peak temperature Tp eak drops almost 
linearly with 7. The slope of this drop is about 1.6 J per 
unit 7. 

In a gaped phase, the low-lying excitation has approx- 
imately the energy dispersion 



£k 



+ ^(k-kof + o((k-k ) 3 ), (2) 



where kg is the wavevector of the excitation minimum, A 
is the energy gap and v the spin velocity. When T <C A, 
it can be shown that xCH nas the form 



X(T) 




(3) 



where A is a T-independent parameter. From the fit of 
the low temperature TMRG results of x(T) with this 
equation, we can estimate the value of A. The result 
of A we obtained is shown in Figure |. The maximum 
energy gap is ~ 2 J/3, located at 7 = 1/3. Our results 
agree with other numerical studies JL7 20 1. 



FIG. 2. The energy gap as a function of 7. 

Figure || shows the temperature dependence of the spe- 
cific heat C(T) for a set of 7. The inset of the figure shows 
the 7 dependence of the peak value of the specific heat, 
Cpcakj and the peak temperature Tp 0ak . Compared with 
^peak' ^peak behaves quite differently It drops with in- 
creasing 7 when 7 < 1/2 and becomes almost a constant 
when 7 > 1/2. Below the peak temperature, C/C pe ak 
shows quite similar behavior for all the curves shown in 
the figure except at very low temperatures. Since there 
is no energy gap at 7 = ±1, C(T) at these two points ap- 
proaches to zero linearly with decreasing T. However, for 
other cases, C(T) decays exponentially at low tempera- 
tures. For the two exact solvable point 7 = ±1, exact 
results are available |2l), the specific heat vanishes lin- 
early at low temperature. However, due to large errors 
at low temperatures, our results do not show this behav- 
ior clearly. Above the peak temperature, C/C pea k drops 
quickly for negative 7. However, when 7 becomes bigger, 
in particular in the incommensurate phase (7 = 2/3 and 
1), C(T) shows a weak and broadened peak above Tp eak . 
It seems that there is a new excitation mode accumulated 
at low energies in the incommensurate state. 

Now let us compare the numerical results with the spin 
susceptibility data Xexp of LiVGe2 0g measured by Millet 
et al. on a powder sample As mentioned in g], two 
extraordinary features appear in Xcxp- Oue is the slow 
drop of Xexp on both sides of the susceptibility peak, and 
the other is the abrupt drop of Xcxp below 22K with a 
sharp upturn below 15K. The first feature, in particu- 
lar the slow drop of Xexp below the peak temperature, 
is reminiscent of a gapless system. The second feature 
of Xcxp is typical of a spin-Peierls system with impuri- 
ties, such as in Zn doped CuGe03 These features 
have led Millet et al. to interpret LiVGe20g as a nearly 
gapless S — 1 spin chain with the spin-Peierls instabil- 
ity. However, whether the abrupt drop of Xexp 

at 22K 

is really due to a spin-Peierls transition is still an open 
question. 
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FIG. 3. The normalized specific heat C/C pe ak as a func- 
tion of T/Tpeak- The upper and lower panels are for 7 larger 
and smaller than zero, respectively. The upper inset shows 
the peak specific heat C pea k and the lower one shows the cor- 
responding temperature Tp eak . J is set to unit. 

The sharp upturn of Xexp at low temperatures indicates 
that the impurity contribution is strong. To see how 
strong the impurity effect is, let us first do a comparison 
without subtracting the impurity contribution in Xcxp- In 
Figure I, the measured susceptibility Xcxp normalized by 
its peak value at about 47 K is compared with the TMRG 
results discussed previously. The disagreement between 
the theoretical and experimental results indicates that 
the impurity effect is too strong to be ignored even at 
high temperatures. 

The susceptibility of dilute magnetic impurities gener- 
ally has a Curie- Weiss form 

c 

Ximp = T + Q' + aT-^ (4) 

where C is proportional to the impurity concentration 
and the square of the effective g-factor of the impurity 
and 8' is a measure for the interaction among impurities. 
The aT~ l term in Ximp is the leading order correction to 
the Curie- Wess term C / (T + 0') due to the finite mag- 
netic field. If there is no interaction between impurities, 
a = (25' 2 + 2S" + l) (g'n B B/k B ) 2 / 10 with g' and S' 
the effective g-factor and spin of impurities. This term is 
not important at high temperatures. But when the tem- 
perature becomes comparable with the level splitting of 
an impurity spin in a magnetic field, this term becomes 
important. It prevents Ximp from being divergent at low 
temperatures, a is typically of order IK 2 when B = IT. 
At very low temperatures the measured susceptibility 



is a sum of Ximp and x(T) given by Eq. (|J), i.e. 

Xcxp(T) = Ximp + A e" A/T . (5) 

By fitting the low temperature experimental data be- 
low 15K with this equation, we find that C" — 
0. 115cm 3 K/mol, 6 1 = U.IK, a = 2.18X 2 , A = 
0.0063cm 3 jmol and A = 36X. These parameters show 
that not only the contribution from impurities to Xexp is 
large as expected, but also the interaction among impu- 
rities is strong at low temperatures. There is no simple 
explanation for such a strong correlation among impuri- 
ties. Clearly this is an important problem which should 
be further investigated both theoretically and experimen- 
tally. 
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FIG. 4. Comparison of the TMRG result (solid line) of the 
spin susceptibility of the S — 1 bilinear-biquadratic model 
with J — 73K and 7=1/6 with the experimental data of 
LiVGe 2 6 §. Xintrin is the experimental data after subtract- 
ing the impurity contribution \imp from Xexp ■ 

By subtracting the impurity contribution from Xexp, 
we obtain the intrinsic susceptibility Xintrin of LiVGe2 06- 
The result of Xintrin together with the raw data Xcxp and 
Ximp is shown in Figure || After the subtraction, the 
abrupt drop of Xcxp at 22K becomes less distinct, but 
the change in the slope is still visible. The most signif- 
icant change of Xintrin compared with Xexp is that the 
peak shifts to a higher temperature and the drop below 
the peak temperature becomes more rapidly. By com- 
paring in detail the normalized Xintrin with the theoreti- 
cal results, we find that Xintrin can be well fitted by the 
numerical curve with 7 = 1/6 (Figure This shows 
that the biquadratic term in model ([j]) does have an im- 
portant contribution to the low energy spin dynamics of 
LiVGe206, in agreement with Millet et al. ||. However, 
the value of 7 which gives the best fit, in particular its 
sign, is different from that suggested in Ref. A de- 
tailed comparison indicates that Xintrin lies between the 
theoretical curves for 7 = 1/4 and 1 /8 in the whole tem- 
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perature region. Thus the uncertainty in the value of j c 
is very small. The result at -f c ~ — 1 suggested in Ref. ||] 
does not fit the experiment data. 

At 7 = 1/6, the peak temperature is Tp Cak = 1.025 J. 
Setting this Tp Cak equal to the peak temperature of 
Xintrin, we find that J ~ 73K. Compared with the gap 
value A = 36K obtained previously, we have A ~ 0.49 J. 
This value of A is rather close to the Haldane gap, 0.54 J, 
of the Hamiltonian ([!]) at 7 = 1/6 (Figure ||) . This sug- 
gests that the low energy spin excitations are gapped and 
the change of the slope at 22K is not due to a spin-Peierls 
transition. 

We have also compared Xintrin with the spin suscep- 
tibility of the S = 1 Heisenberg model with uniaxial 
single-ion anisotropy but without the biquadratic term 
@, namely the model fl" = J£< [S< • S i+1 + D £\ Sf z ] . 
However, in the parameter region which might be phys- 
ically relevant, —1/2 < D < 1/2, we find that none of 
the numerical curves fits Xintrin in the whole tempera- 
ture range. This shows that the uniaxial anisotropy in 
LiVGc20g is indeed very small, in agreement with the 
analysis of Millet et al. . 

The above analysis confirms the importance of the bi- 
quadratic exchange interaction in LlVT^Og. On the 
other hand, it also raises some new questions. In the 
argument given by Millet et al., the biquadratic term 
comes from in fourth order since at second order the an- 
tiferromagnetic and ferromagnetic terms partially can- 
cel. However, the coefficient of this biquadratic term is 
negative (i.e. 7 < 0) according to their calculation, in 
contrast with the result we obtain. To resolve this dis- 
agreement, further investigation into the electronic struc- 
ture of LiVGe2C>6 is needed. More detailed measure- 
ments with high quality single crystals would also help 
clarify the impurity effect as well as the nature of the 
anomaly at 22K in this material. In a S — 1 Heisenberg 
chain, the localized non-magnetic impurity may induce 
mid-gap states within the Haldane gap [^4|^7| . A better 
understanding of the physical properties of these mid-gap 
states would also be helpful for further understanding the 
thermodynamics of LiVGe2 06 at low temperatures. 

In summary, the thermodynamic properties of the bi- 
linear and biquadratic Heisenberg model have been stud- 
ied and compared with the experiments. The measured 
susceptibility data of LiVGe 2 Og, after subtracting the 
impurity contribution, can be quantitatively explained 
by the model (Q) with 7 = 1/6. 
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